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Abstract. This study introduces a modified quadratic Lorenz 
attractor. The properties of this new chaotic system are analysed 
and discussed in detail, by determining the equilibria points, the 
eigenvalues of the Jacobian, and the Lyapunov exponents. The nu¬ 
merical simulations, the time series analysis, and the projections 
to the xy-plane, x^-plane, and r/ 2 -plane are conducted to high¬ 
light the chaotic behaviour. The multiplicative form of the new 
system is also presented and the simulations are conducted using 
multiplicative Runge-Kutta methods. 


1. Introduction 

Dynamical systems are mathematical models describing the evolu¬ 
tion of systems in terms of equation of motion with sensitive initial 
values. The application areas of dynamical systems can be found e.g. 
in many disciplines in physics [Ml DM, population dynamics in biol¬ 
ogy [32], and chemical kinetics in chemistry [29i|25]. Dynamic systems 
theory has also found its way into subjects outside the fields of Math¬ 
ematics and natural sciences, exemplarily we would like to mention 
mathematical economy and finance mmm- Poincare made a kick 
start to the subject of dynamical systems in his pioneering work in 
1890 [2T] . Later, in the 1920’s, Fatou [8] and Julia [ID] introduced the 
dynamics of complex analytic maps. On the one hand, studies like 
Birkhoff [3], Kolmogorov [T3J Mi IJ2J, Cartwright and Littlewood [3] 
find their origin in physics, like e.g. the three body problem in as¬ 
tronomy, whereas on the other hand Stephan Smale provided a purely 
mathematics motivated approach [271 EH]. E. N. Lorenz [16] observed 
that very simple differential equations become chaotic under certain 
circumstances. The system proposed by Lorenz shows a very complex 
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dynamical behaviour and displays the well-known two-scroll butterfly- 
shape. The dynamic equations of the Lorenz system are given as 

(1) d i = 

( 2 ) ^ = x(p^z)-y, 

dz 

(?) Jt = xy-fr, 

where the parameters a, p , and f3 are assumed to be positive. Lorenz 
used exemplarily the values a = 10, /3 = 8/3 and p = 28 to demon¬ 
strate the systems chaotic behaviour. The study on ’’The equation for 
continuous chaos” by Rossler [23] can be seen as another landmark in 
the discussion of 3D dynamic systems. The Rossler system and the 
Burke Shaw system [25] have the property of two unstable saddle foci 
in common. Other chaotic systems exhibiting a similar simple struc¬ 
ture as the Lorenz system, without being topologically equivalent, are 
proposed by Chen [5| and Lii and Chen m ca discusses the tran¬ 
sition between Lorenz and Chen attractors. Furthermore, Yang and 
Chen introduced another chaotic system with three fixed points: one 
saddle and two stable fixed points. Yang et al. [30] and Pehlivan et 
al. [19] introduced and analysed chaotic systems similar to the Lorenz, 
Chen, and Yang-Chen systems, with two different fixed points, i.e. two 
stable node-foci. Chaotic systems have found their ways also into many 
applications in engineering, such as electronic circuits E3-EB- Modern 
pacemakers actually rely nowadays on these chaotic circuits. 

Dynamical systems have also been discussed in the framework of var¬ 
ious Non-Newtonian Calcului, like fractional calculus, geometric mul¬ 
tiplicative calculus, and bigeometric multiplicative calculus. Jun Guo 
Lu transformed the Lii system into fractional calculus [18], investigat¬ 
ing chaotic behaviour of fractional-order of the Lii system numerically. 
Aniszewska applied multiplicative calculus to the Rossler system and 
showed chaotic behaviour of multiplicative Rossler [I]. As an applica¬ 
tion of the bigeometric multiplicative Runge-Kutta method, the bige¬ 
ometric multiplicative Rossler System was solved in [23]. 

This study introduces a new chaotic attractor, found by modification 
of the Lorenz system by a quadratic term. Detailed numerical and 
theoretical analysis reveals that the proposed system shows chaotic 
behaviour and the property of a two-scroll attractor like the Lorenz 
attractor. 

Section [2] introduces the modified quadratic Lorenz attractor, dis¬ 
cusses and analyses its properties by determining the equilibria points, 
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and Lyapunov exponents theoretically as well as numerically. In Sec¬ 
tion 3, the modified quadratic Lorenz attractor is translated into geo¬ 
metric and bigeometric calculus, and the solutions of the the system are 
obtained using the corresponding multiplicative Runge-Kutta methods 
[22 J and |23j . Finally, this paper closes with the summary of the ob¬ 
tained results. 


2. Design of a new Chaotic System 


This paper presents a new Chaotic system derived from the Lorenz 
system. The system is generated by the following simple three-dimensional 


system: 



(4) 

dx 

= <y(yz - x) 

dt 

(5) 

dy_ 

dt 

= px — xz 

(6) 

dz 

dt 

= (xy f - (3z 


where x, y, and z are variables and er, p, and (3 are real parame¬ 
ters. In the new proposed chaotic system, all the equations have some 
differences compared to the original Lorenz system. In order to see 
the differences between the two systems we can compare the equations 
0-0 and @)-@ one by one. Evidently, in the Lorenz system equa¬ 
tion ^l| is linear, whereas equation Q is non-linear. Furthermore, 
equations p]) and ([5]) are both non-linear, where the y-dependence of 
equation pF is cancelled. The most significant difference can be ob¬ 
served comparing equations ([3]) and (j6j) . In Q the term xy is squared 
compared to equation (|3]) . 


2.1. System Description. The initial values and the parameters of 
the system are chosen as (1,1,1) and a = 12, p = 8 with varying (3. 
Then it has been observed that the behavior of the new chaotic system 
changes depending on the different values of j3 as shown in Figure [lj 
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(c) (3 = 2 (D)/3 = 4 



Figure 1. Simulation of the system when a = 12, p = 8 
and various parameter values 

As a result of the calculations, as it can be seen in Figure [lj it has 
been observed that the new system is chaotic for the parameters 

(7) a — 12, p = 8, and (5 = 4. 

The analysis of the new chaotic system will be done according to those 
parameters. 

2.2. System Analysis. The first step to analyze a chaotic system is 
to find the equilibrium points. In order to determine the equilibrium 
points of the proposed system 0-& we need to solve the system 

a(yz - x ) = 
px — xz = 

(xy) 2 - pz = 


(8) 


0 

0 

0 
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Thus, the solution of the system ^ with respect to x, y, z give the 
equilibria points as: 


O = ( 0 , 0 , 0 ) 



For the parameters chosen in (|T]) , the numerical values of the equi¬ 
libria points are 


(9) 

( 10 ) 
( 11 ) 


0 = ( 0 , 0 , 0 ), 

E + = (6.73,0.84,8), 

E~ = (-6.73,-0.84,8). 


In order to decide on the stability of the new proposed system, the 
eigenvalues of the Jacobian matrix must be analyzed. The Jacobian 
matrix for this system (|4])-([6]) can be easily obtained as 




—a 

<JZ 

ay 

(12) 

J = 

p-z 

0 

—x 



2 xy 2 

2 x 2 y 

-P 


The expressions for the eigenvalues of the Jacobian matrix (12) are 
very long and complicated. As we are only interested in the numerical 
values of the eigenvalues at the equilibria points (|9])-([TTj) for the given 
parameters ([7]), the eigenvalues are stated in the table below: 


Equilibrium 

Point 

Ai 

a 2 

A 3 

O 

-12 

-4 

0 

E + 

2.65 + 23.87* 

2.65 - 23.87* 

-21.3 

E~ 

2.65 + 23.87* 

2.65 - 23.87* 

-21.3 


TABLE 1 . Eigenvalues of the Jacobian at the equilibrium 
points 


As the eigenvalues Ai and A 2 for the equilibrium point O are both 
negative, the system is unstable at this equilibrium point. The eigen¬ 
values corresponding to the equilibrium point E~ will be the same with 
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the eigenvalues of E + , because of the quadratic nature of the system. 
Since A 3 is a negative real number and Ai and A 2 are two complex con¬ 
jugate eigenvalues with positive real parts, equilibrium points E + and 
E~ are unstable according to [9]. 

2.3. Symmetry and Dissipativity. The System (|4])-(|6]) has a nat¬ 
ural symmetry and is invariant under the coordinate transformation 
(x,y,z) —> (— x,—y,z ) which persists for all values of the system pa¬ 
rameters. So, system Q -© has rotation symmetry about the z — axis. 

Let, /1 = , fi = J and / 3 = § in the system Q-Q. Then we 

get for the vector held 


(x,y,z) T = (/i,/ 2 ,/ 3 ) T 

Consequently the divergence of the vector held V yields to: 

V7 r ■ -\t dfi df 2 df 3 

Note that / = — (cr + /3) = —16 is a negative value, so the system is 
a dissipative system and an exponential rate is: 


(14) 




V(t) = Voe* = V 0 e 


- 16 1 


From (14), it can be seen that a volume element Vo is contracted by 
the how into a volume element Voe~ 16t at the time t . 


2.4. Lyapunov Exponent and Fractional Dimension. The Lya¬ 
punov exponents generally refer to the average exponential rates of 
divergence or convergence of nearby trajectories in the phase space. 
The important part is that if there is at least one positive Lyapunov 
exponent, the system can be defined to be chaotic. According to the de¬ 
tailed numerical and theoretical analysis, the Lyapunov exponents are 
found to be l\ = 5.4162, l 2 = 2.1912, and Z 3 = —19.2269. Therefore, 
the Lyapunov dimension of this system is: 


(15) 




D L =i + AA =2+ h+h 


l(j+l 


2.3957 


This result is consistent with the findings in jUj, i.e. that the Lyapunov 
dimension is in the range 2-3. Equation © means that the system 
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is a dissipative system, and the Lyapunov dimensions of the sys¬ 
tem are fractional. Having a strange attractor and positive Lyapunov 
exponent, it is obvious that the system is a 3D chaotic system. 


Dynamics of Lyapunov exponents 



FIGURE 2. Plot of Lyapunov exponents 


2.5. Numerical Simulations. Time series analysis of the system Q- 
(|6]) according to x(t),y(t), z(t) axes are listed in the Figure [3j For the 
solutions of the system the Runge-Kutta method was employed. 
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Obviously the time series show that the functions x(t), y(t), and z(t) 
are not periodic, indicating that the system is chaotic. 

The projections of the system (|4])-(j6]), on the various axes is given in 
the Figure |4j 
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(C) 

Figure 4. Projection of System 0-0 on the x-y 
plane, x-z plane, y-z plane respectively 


The chaotic attractors are displayed in the Figure [4j It appears that 
the new attractor exhibits an interesting complex chaotic dynamics 
behavior. 

3. Geometric and Bigeometric Sense of New System 

As an application of the dynamical systems, the newly formed chaotic 
system can also be written in the sense of Geometric and the Bigeo¬ 
metric calculus. 

3.1. Geometric and the Bigeometric Chaotic Systems. The Geo¬ 
metric Chaotic system corresponding to system 0-0 , can be written 
as 


d*x 

dt 

d*y_ 

dt 

d*z 

dt 


cr(yz-x) 

e x 


px — xz 

e v 


(xy) 2 -l3z 

e 


(16) 
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On the other hand, the Bigeometric Chaotic system corresponding 
to system Q-Q will be 


(17) 


OC X t(a(yz-x)) 

—— = e x 

dt 

d n y t(px-xz) 

—— — e y 
dt 

d n Z t((xy) 2 -pz) 

—— = e 
dt 


In order to solve the system (16), one can use the Geometric Runge- 
Kutta method defined in [(22] . while the system 0 can be solved by 
using the Bigeometric Runge-Kutta method which was defined in [23]. 
By choosing the same values for the parameters, such as cr = 12, p = 8 
with varying f3 and the initial values as (0,0,0), we can see that the 


solutions of the Geometric chaotic system (16) and the Bigeometric 


chaotic system (17) will be similar to the ones that we get from the 
solutions of the chaotic system 0-@. The following figure shows 
the chaotic behaviour of the Geometric and the Bigeometric dynamic 
systems when the parameters are chosen as a = 12, p = 8 and (3 — 4. 


z 



Figure 5. 0 = 4 
4. Conclusion 

The proposed modified quadratic Lorenz attractor was analysed the¬ 
oretically and numerically showing that this system is chaotic. Follow¬ 
ing the standardized analysis method for chaotic systems, we deter¬ 
mined first the equilibria points and the eigenvalues of the Jacobian 
matrix at these equilibria points to get a first indication about the 
stability of the proposed system. As the eigenvalues are either non¬ 
positive real numbers, or complex numbers with positive real parts, 
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we can conclude that this system is not stable at the equilibria points. 
Moreover, we could identify that the proposed system has a rotational 
symmetry about the 2 -axis, and shows dissipative behaviour contract¬ 
ing the volume element Vo to Voe~ 16 . The Lyapunov exponents yield to 
l\ = 5.4162 ,12 = 2.1912 and Q = —19.2269, showing the chaotic nature 
of the system. The fractional dimension of the system Dl — 2.3957 has 
also been given. Overall, the analysis shows that this is a new chaotic 
system with two scrolls. 
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